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Abstract 

An ultra-stable optical clock based on coherent population trapping effect of alkaline-earth ions, 
such as Ca^, Sr^, Ba^, is analyzed here. The proposed transitions use the odd isotopes, so that the 
frequency shift is insensitive to the applied magnetic field. To enhance the signal, a large number 
of ions are trapped in a linear Paul trap, and laser cooled to crystallize such that the ions are in 
the Lamb-Dicke regime to avoid first order Doppler shift. Other relevant frequency shifts are also 
analyzed. Further techniques to improve the stability of the clock signal are discussed. 
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Optical clocks, or atomic clocks using optical transitions, have much higher line Q com- 
paring to microwave atomic clocks. As a consequence, optical clocks potentially can have 
better stability and accuracy than today's microwave atomic clocks. Many groups are work- 
ing on optical clocks with the goal of realizing improved frequency standards and testing 
fundamental physics. Commonly, there are two ways of making an optical clock. One is 
based on a single, tightly trapped ion with very small systematics uncertainties. The dis- 
advantage of this method is the intrinsically smaller signal strength. Consequently, a long 
sampling time will be required to realize the potential of this type of optical clocks. This 
problem can be alleviated by confining many neutral atoms in an optical lattice trap with 
much improved signal strength [ij]. In this case, the relevant atomic levels shift equally (via 
AC Stark shift) in the optical lattice, at the so-called "magic wavelength." Both approaches, 
so far, operate in pulsed interrogation mode. Pulsed interrogation, such as "shelving," has a 
low duty cycle that relies heavily on the local oscillator's free-running stability. In addition, 
there could be complicated aliasing problems such as the Dick effect j^. 

In this Letter, we analyze a third approach, namely, an optical clock that operates contin- 
uously using the coherent optical trapping (CPT) transitions of a large, laser-cooled, trapped 
ion crystal. We consider a cloud of ions trapped in a large, linear ion trap. In this case, the 
ions will be cooled to sufficiently low temperature such that an ion crystal is formed. At 
such low temperatures (~ 1 mK), the ions' secular motion amplitude is less than half of the 
laser wavelength, and we are in the so-called Lamb-Dicke regime I2I , resulting in a vanishing 
first-order Doppler shift. This situation is identical to the optical lattice, where an artificial 
optical potential is applied to confine the atoms. In the case analyzed here, simple laser 
cooling will crystallize and confine the ions, without the additional lattice laser. Under such 
conditions, it is possible to construct an optical clock on the sharp dark lines of CPT with 
vanishingly small systematics. 

CPT, a phenomenon that is accomplished by two coherent lasers interacting in a A scheme 

nn 

[^,15|, is generally observed in alkali atoms. There are many applications with CPT, including 
microwave atomic frequency standards 6] and magnetometry Q, Bl Narrow linewidths 
below 50 Hz have been realized in buffer-gas-filled vapor cells |9|, [lO|. This corresponds 
to relative clock frequency stability as low as 1.3 x 10^^^ for 1 s integration time. Even 
better stability can be expected with an optical CPT. In this Letter, we propose an optical 
frequency standard based on optical CPT, as shown in Fig. 1(a). The proposed scheme is 
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applicable to alkaline-earth ions, such as Ca"^, Sr+, and Ba"^. But as an example, we will 
focus on Ba"*". In the proposed scheme, two phase-locked lasers at 493 nm (green laser) and 
650 nm (red laser) will interact with crystallized Barium ions, and spectrally narrow light 
at the exact clock frequency can be generated through frequency mixing as the frequency 
standard. 

More specifically, odd isotope ^^^Ba"*" will be chosen, and the three energy levels involved 
are |65'i/2, F — 2), |6Pi/2, F — l), and |5£)3/2, F — O), as shown in Fig. 1(b). It is possible 
then to use circular polarized light to estabhsh coherence between |65'i/2,-F = 2,mF = O), 
and I5D3/2, = 0, mp — O). The transition between mp — ^ mp = is chosen so that 
they are independent of the applied magnetic field in first order, a desirable property in 
frequency standards applications. 

To understand the full dynamics of the proposed system under magnetic field, a twenty- 
level density matrix would have to be used. But the basic phenomenon that we are interested 
in can be understood by using a three-level model. Under the rotating wave approximation, 
the time evolution of the three- level atomic system coupled to the two hght fields is governed 
by the Liouville's equation, 

| = (1) 
where p is the atomic density matrix, and we have 

Tr[p] - pn + P22 + P33 = 1- (2) 



The system Hamiltonian is defined as 



^ Ag Q12/2 ^ 



(3) 



Q12/2 Q23/2 
y Q23/2 A, J 

where = Ug — Ui2 and A^ = a;^ — a;23 are the detunings between the green (|1) |2)) and 
red (|2) |3)) laser frequencies, cUg and cUr, and the atomic resonance frequencies, UJ12 and 
UJ23, respectively. Q12 and ^,23 are the Rabi frequencies associated with the couplings of the 
green laser and the red laser to atomic transitions |1) — > |2) and |2) — > |3), respectively. The 
term L(p) includes atomic relaxation terms and laser linewidths effect, and can be written 
as 

Up) = -I J2(^nCnP + pClCn - 2CnpCl). (4) 



2 

n 



Since 5D3/2 state has a very long lifetime (80 s), we can neglect its decay to 65'i/2 state. 
The decay from 6P1/2 state to 66*1/2 state is described by 

(5) 



C21 = vr2i|i)(2|, 

where r2i = 27r x 15.1 MHz. Similarly, the decay from 6P1/2 state to 5D3/2 state is described 
by 

C23 = v^|3)(2|, (6) 
where r23 = 2tt x 5.3 MHz. The effect of laser linewidths can be introduced by 

= v^|l)(l|, (7) 

(8) 



Cr = v^|3) (3| 



Here, Tg and Tr denote the laser linewidths of the green and red lasers, respectively. Sub- 
stituting Eq. ([5]), Eq. ([6]), Eq. (JTj), and Eq. (IHl) to Eq. (jlj), we obtain the Liouvilhan matrix 

Hp), 

^ r2iP22 -(r2i + r23 + 2r3)pi2/2 -(rg + r,)pi3 ^ 

L{p)= -(121 + r23 + 2r,)p2i/2 -(r2i + r23)p22 -(r2i + r23 + 2r,)p23/2 

-{Tg + Tr)p31 -(121 + r23 + 2rr)p32/2 r23P22 



V 



(9) 

Under stationary condition where p = 0, Eq. ([T]) can be solved to obtain the resonance 
fluorescence excitation spectrum under different detunings. When = A^, a CPT dark 
line will appear, as shown in 11, Q, [l^. The linewidth of the dark line depends on laser 
linewidths and laser power. If we neglect the effect of laser linewidths, then the CPT 
linewidth at optical resonance is 



Au = -^==\ (1 + 3S + -S*){1 + -S*) 



(10) 



Here, 7 is the decay rate of the 5D3/2 state, F = r2i + r23, and S and S* are two saturation 
parameters, 

^12^23 ^ 



S 



s* 



^12 "I" ^23 
Ql2 + ^23 



4r2 



(11) 

(12) 
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At zero laser powers, the CPT linewidth is hmited by the decay rate of the 5D3/2 state at 
2 mHz. More reahstic value of ~ 10 Hz can be achieved with reasonable laser powers. As 
an example, using ni2/27r = Q23/27r = 10 kHz, the calculated linewidth is 5 Hz, as shown 
in Fig. 2(a). Experimentally, the peak of this dip can be locked by dithering the driving 
frequency of an AOM that is used to scan the repump laser frequency. Assuming 10^ ions 
with 1% collection efficiency, we can achieve a PMT counting rate of around 10 kHz, from 
which we can estimate that the proposed frequency standard has stability of better than 
10"^^ in 1 s in the shot-noise-limited photon counting situation. 

At this moment, it is worthwhile to estimate all systematic frequency shifts of the pro- 
posed clock transitions. In second-order perturbation theory, the light shift of state \k,m) 
under light beam of frequency a;, electric field strength F, and polarization e can be calcu- 
lated as 

k' ,m' ,±w 

Here, k stands for atomic quantum number nJF, and is the unperturbed energy of level 
k. The sums extend over continuum states. In our case, for |65'i/2, F — 7^ state, the main 
light shift contributions come from the near-resonant transition light of 493 nm connecting 
|65i/2, F = 2> to |6Pi/2, F = 1> and |6Pi/2, F = 2> states; for |5F'3/2, F = 0> state, the main 
light shift contributions come from the near-resonant transition light of 650 nm connecting 
|5Z^3/2, F = 0) to |6Pi/2, F = 1) state. As a consequence, we can estimate the light shift of 
the |6S'i/2, F = 2) state under right-circular polarized light as 

47rl]2^ l(65i/2,F = 2,m = 0|/il|6Pi/2,F' = l,m' = l>|' 



3eoA3r2i ' hAg 

I (6^1/2, F = 2,m = 0|/xl|6Pi/2,F^ = 2,m' = l)f 

h{Ag-6p) J' ^'^^ 

where Sp/2'K = 1488 MHz is the 6P1/2 state energy hyperfine splitting. Here, we have used 
the fact that 

Ig = ^ceoFj 
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The two matrix elements in Eq. (fT^ can be calculated as 

|(6^i/2,F = 2,m = 0|^}|6Pi/2,F' = l,m' = l)| = -i=/i„ (16) 
|(65i/2,F = 2,m = 0|/i}|6Pi/2,F' = 2,m' = l)| = (17) 

where the reduced matrix element jjLg = |(6S'i/2| I6P1/2) | = 8.43 Debye. Similarly, the 
light shift of the \5D^/2,F = O) state under right-circular polarized light is 

Anill, |(5D3/2,i^ = 0,m = 0|/il|6Pi/2,P' = l,m' = l)|' 

where 

|(5Z}3/2,P = 0,m = 0|/xi|6Pi/2,P' = l,m' = 1)1 = -i=/i„ (19) 

with the reduced matrix element jjr = | (5-03/21 lAi"*^! I6P1/2) | = 7.64 Debye. By adjusting the 
relative intensities of the green and red lasers, it is possible to suppress the light shift to a 
negligible level. For example, using the parameters of Fig. 2(b), with f2i2/27r = 14 kHz and 
f223/27r = 10 kHz, the total light shift can be calculated as 

Auls = Au5D - Az/6s = -4.8 mHz. (20) 

This value can be further decreased if we choose greater detunings. With the parameters of 
Fig. 2(b), the transition linewidth is increased, but the signal strength is similarly increased, 
so the resulting clock stability is still effectively the same. 

For the proposed clock transition, the first-order Zeeman shift is zero, and the quadrupole 
DC-Stark shift due to stray electric field gradients is also zero in first-order. The accuracy of 
the proposed frequency standard is primarily limited by second-order Doppler shift, second- 
order Zeeman shift, blackbody Stark shift, and quadrupole AC-Stark shift. 

The second-order Doppler shift can be estimated as 

Assuming a secular frequency uJs/^tt = 0.1 MHz, and the mean occupation number (n) = 
100, the second-order Doppler shift can be estimated to be —1.65 x 10"^^. Using Breit-Rabi 
formula p^], the second-order Zeeman shift can be calculated as 



With an applied 1 mG magnetic field and a 0.1 mG stability, the second-order Zeeman shift 
is —7.55 X 10^^^ with a relative uncertainty of 7.55 x 10"^^. The blackbody Stark shift can 
be calculated as 



^^BB _ ^BB,5D — ^BBfiS 
V V 

n 

Here the time averaged electric field is pJi] 



6.37 X 10-^^(r/cm)-2 (e^) . (23) 



E^) = (8.319 V^/cm)^ (^^^ . (24) 

At 300 K room temperature, we can estimate the blackbody Stark shift to be —5.1 x 10^^^. 
With a temperature stability of 1 °C, we estimate an uncertainty of the shift at 6.8 x 10^^'^. 
The 6S1/2 quadrupole AG-Stark shift is zero, while 5D3/2 quadrupole AG-Stark shift is 
estimated to be 0.2 mHz if we assume 1 mm trap size and 200 V rf voltage, so the final 
quadrupole AG-Stark shift is 

^^Stark ^Stark,5D ^ -j^g—ig f25) 

Hence, we can see that the dominant systematics is due to blackbody Stark shift. How- 
ever, its influence on clock stability can be eliminated by better control of the surrounding 
temperature variation. 

In order to realize the proposed frequency standard, there are several practical issues that 
need to be dealt with. First, the GPT process connect |6S'i/2,-F = 2), |6Pi/2,-F = l), and 
|5D3/2,-F = 0) states, but |6Pi/2,-F = l) can also decay to 1651/2,-^ = l), |5-D3/2,P = l), 
and |5D3/2,-F = 2) states. This problem can be solved by using another set of auxiliary 
cooling and repumping lasers to close the pumping cycle. Second, since we use circular 
polarized light, atoms are optically pumped to a Zeeman sublevel of |6S'i/2,-F = 2) that 
is not involved in the GPT process. This effect has consequences on the observed signal 
amplitude. To fully account for this effect, an effective four level model would have to be 

n 

use d | 16l| . This effect of optical pumping can be reduced by using a double lambda scheme 
18| . Finally, in order for the proposed GPT scheme to work, the laser linewidths need to 
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be stabilized to the Hz level. This is readily achievable with current development of ultra- 



stable high finesse reference cavities 19|, l20|, l2l|, |22]. Also, precise relative phase coherence 



between the green and red lasers are required. This can be achieved by locking a frequency 
comb to one of the laser, serving as the master laser, and phase-lock another laser to the 
frequency comb with a tunable frequency offset. 



To further improve the stabihty of the frequency standard, we need to further decrease 
the CPT hnewidth and increase the signal strength. This is difficult to do in continuous 
CPT due to optical power saturation. An alternative method to overcome this limit is to 
generate optical CPT pulses to observe Ramsey fringes 23|, |2J]. With this method, the 
signal strength can be increased greatly, while the Ramsey fringe widths are dependent only 
on the transit time T between the interaction zones as 1/2T. This time domain Ramsey 
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27|, and we believe it will also 



method has been demonstrated in microwave CPT 
work for optical CPT. 

In conclusion, we propose and analyze in details an ultra-stable optical frequency stan- 
dard based on optical CPT. The proposed frequency standard can work in the continuous- 
wave mode and has a frequency stability at the 10~^^ level, with very favorable systematics 
uncertainties. Even better stability can be achieved with Ramsey CPT methods. 

ZHL acknowledges fellowship support from Alexander von Humboldt Foundation. 
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(a) (b) 

FIG. 1: (Color online) (a) Energy level structure and optical couplings of the three-level atomic 
system for making a CPT frequency standard, (b) Energy level structure of odd isotope Barium 
ion. 
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FIG. 2: Resonance fluorescence excitation spectrum with different red laser detuning. Assuming 
10^ ions and 1% collection efficiency, (a) Ag/27r = -20 MHz, rii2/27r = J^23/27r = 10 kHz, 



= 0. (b) Ag/2'!T 



-20 MHz, J^i2/27r = 14 kHz, n23/27r = 10 kHz, Fg = 



0. 
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